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We develop an analytical theory of polarization of a vertically aligned array of carbon nanotubes 
(NTs) in external electric field. Such arrays are commonly utilized in field-emission devices, due to 
the known electrostatic effect of strong field enhancement near the tip of an individual NT. Small 
ratio of the NT radius to the separation between neighboring NTs allows to obtain asymptotically 
exact solution for the distribution of induced charge density along the NT axes. For a regular array, 
this solution allows to trace the suppression of the field penetration with increasing the density of 
NTs in the array. We demonstrate that for a random array, fiuctuations in the NT density terminate 
the applicability of our result at distances from the NT tips much larger than the field penetration 
depth, where the induced charge density is already exponentially small. Our prime conclusion is 
that, due to collective screening of the external field by the array, the field-emission current decreases 
drastically for dense arrays compared to an individual NT. We argue that the reason why the strong 
field emission, described by the Fowler-Nordheim law and observed in realistic arrays, is the strong 
dispersion in heights of the constituting NTs. 

PACS numbers: 73.40.Gk, 79.70.+q, 79.60.Ht, 81.07.De, 85.35.Kt 



I. INTRODUCTION 

First report of field emission from carbon NTs 
had appeared a decade agoi. It was followed by a 
demonstration^ that arrays of NTs can be patterned into 
emitting and non-emitting regions. Since then, the field 
emission properties of carbon NTs command a steady 
interest from researchers worldwide. The uniqueness 
of these properties originates from geometry of a NT. 
Namely, due to a small NT radius, r, the electric field 
applied between the substrate (cathode), on which the 
NTs are grown, and the anode is enhanced by a large fac- 
tor /3 3> 1 near the nanotube tip. Such an enhancement 
translates into high probability of electron tunnelling to- 
ward the anode, leading to desirable low turn-on voltage 
for field emission. This property, combined with high 
emission current density, made possible a successful fab- 
rication of the row-column matrix-addressable flat panel 
display based on carbon NTo^i'^i^i'^/i^i^d°d^ Currently, 
fiat panel displays constitute one of the most prominent 
applications of NTs^^. Geometrical characteristics of in- 
dividual NTs utilized in the first displajii were highly 
dispersed. Further advances in fabricationi^ allowed to 
achieve excellent vertical alignment and high homogene- 
ity in the lengths and radii of nTe^'*'!''^-!'^-!'^. This sug- 
gests that NTs on the cathode of a display can, in the 
first approximation, be viewed as constituting a regular 
array of identical NTs. Such an array is schematically 
illustrated in Fig. 1. 

On the theoretical side, the focus of the previous 
studieoi«d9,20,2i,22,23,24^25,26,27 ^^^A cmissiou from nan- 
otubes was the effect of band structure and tip geometry 
of an individual NT on the emission probability. These 
studies left out the fact that all NTs are coupled to each 
other electrostatically. 

The main point of the present paper is that, for a reg- 
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FIG. 1: Schematic illustration of the forrest of NTs of a 
height, h, grown on the substrate, z — 0. The average dis- 
tance between neighboring NTs is d. Shaded regions at the 
NT tips illustrate the charges, induced in the NTs by the ex- 
ternal field, F. For a dense forrest, the typical penetration 
depth, a, exceeds d. Vertical arrows illustrate electric field. 
Find, created by induced charges. Regions of lower NT density 
correspond to higher field-enhancement-factor, /3a = Fi„d/F. 
In these regions, external field penetrates deeper into the for- 
rest. Fluctuations in the induced charge density, due to the 
randomness in the NT positions, exceed the average density 
at depth {h — z) > a'^ /d'^. 

ular NT array mutual electrostatic coupling of NTs in 
the array has a dramatic effect on the field emission, es- 
pecially in dense arrays. By dense we mean the arrays 
in which the separation, d, between neighboring NTs is 
much smaller than the NT height, h. The situation d <^h 
is quite common in realistic field-emission setups, see, 
e.g., R,efs. IT^HT^ItR To substantiate this point, consider 
first two parallel NTs separated by d <C in the external 
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electric field, F, directed along their axes, see Fig. 1. It 
might seem, that, if d exceeds the tunnelling length for 
field emission from each of the NTs, then both NTs emit 
independently. This, however, is not the case. The rea- 
son is that the enhancement of the electric field near the 
tip of each NT is governed by the charge density, induced 
by the external field. For two parallel NTs, the induced 
charge density per NT is approximately two times smaller 
than for an isolated NTi2& As a result, the field enhance- 
ment near the tip of each NT becomes smaller due to the 
presence of the neighbor. 

Compared to a pair of neighboring NTs, the suppres- 
sion of the field enhancement becomes much more pro- 
nounced in the NT array. On the qualitative level, this 
conclusion follows from the fact that each NT in the ar- 
ray interacts with ~ (h/d)'^ ^ 1 neighbors. As we will 
demonstrate below, external field simply does not pene- 
trate into the sufficiently dense array. On the intuitive 
level, this trend has been previously understoodiiS422i2£ 
In particular, in Ref. numerical simulations illustrat- 
ing the suppression of the field enhancement for three 
parallel NTs with decreasing distance between them were 
reported. However, full understanding of this suppression 
in the array with arbitrary ratio d/h requires an analyt- 
ical theory. Such a theory is developed in the present 
paper. In particular we demonstrate that 

i) penetration of the external field into the array is de- 
scribed by a simple function, sinh(z/a), where the pene- 
tration depth, a, is much smaller than h for a dense array, 
but still much bigger than d; 

ii) distributions of induced charge density in a regular 
and completely random NT arrays are approximately the 
same] 

Hi) with regard to the field emission, the enhancement 
of external field for the array, as compared to the indi- 
vidual NT, is suppressed by the factor k, [h/a). 

The reason why the electrostatic problem, in which 
the variables cannot be separated, allows an asymptotic 
analytical solution is the presence of large parameters 
{h/r) ^ 1 and {d/r) ^ 1. As it was demonstrated in 
Ref. Isil, for a single NT, the relation h ^ r allows one 
to find analytically the distribution of induced charges in 
external field. Here the approach of Ref. ^Jjis generalized 
to the NT array. 

The paper is organized as follows. In Sect. II we re- 
view the Thomas-Fermi description?^ of polarization of 
a single NT in external field. In Sect. HI we generalize 
the Thomas-Fermi equation to the NT array. In Sec- 
tions IV-V we analyze this equation for a regular array, 
and find its asymptotic [in the parameter (d/r) 3> 1] 
solution. Robustness of this solution with respect to ran- 
domness in the NT positions is demonstrated in Sect. VI. 
In Sect. VII we apply to the obtained solution for distri- 
bution of the induced charge density to calculate the field 
emission current from the array. Relation of our theory 
to experiment is addressed in Sect. VIII. 



II. SINGLE NT 

Denote with p{z) the linear charge density on the NT 
surface at a distance z < h from the substrate. The 
Thomas-Fermi equation for p(z) reada^i 

1 1 

eFz=-p{z) + - dz' So{z,z')p{z'), (1) 
9 e 7o 



where the kernel, Sq{z,z'), is defined as 

Sa{z, z') = ^{z - z') - <^{z + z'). 



(2) 



Here the function $(z) is the interaction potential be- 
tween two points on the NT surface, separated vertically 
by z. It represents an azimuthal average of the Coulomb 
potential 
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(3) 



where r is the NT radius. The second term in Eq. (0) 
accounts for the image charges. The meaning of Eq. 
is the following. The Ihs is the bare potential. The po- 
tential, acting on a given electron at the NT surface rep- 
resents the sum of this potential and of the potential, 
created by the induced charges. This resulting potential 
defines the local value of the Fermi energy, which, in turn, 
fixes the local value of the Fermi momentum. This Fermi 
momentum, on the other hand, is linearly proportional to 
the charge density in one dimension. This is a standard 
reasoning behind the Thomas-Fermi description. Within 
this description Eq. is nothing but the condition, that 
the electrochemical potential remains constant along the 
NT. 

The prime simplification that allows analytical solution 
of Eq. is that, with logarithmic accuracy, p{z') in the 
integrand in the rhs of Eq. ^ can be replaced by p{z). 
Upon this replacement, we have 

If'', , . h^ ( l-z/h\ , , 

-j^ dz 5o(z,z)«ln--ln--Kln(^^^j, (4) 

where we assumed that {h — z) 3> r. With the same 
logarithmic accuracy, for z ^ r, the rhs of Eq. Q can be 
replaced by 2£;i, where = hi{h/r). Then we obtain 
the following analytical solution for the induced charge 
densitjiSi 



p{z) 



gFz 



1 -I- 2gCh 



(5) 



where we have introduced a dimensionless interaction pa- 
rameter g — 2Ne'^ /nhvQ. The above result for p{z) is ap- 
proximate, in the sense, that the numerical factor in the 
argument of a large logarithm, iZ^ 3> 1, is not specified. 
Equation Q represents the solution of Eq. (Q, which 
satisfies the obvious condition p(0) = 0. An improved 
analytical description yielding the result coinciding with 
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Eq. ISI) in the limit of large Ch was recently reported in 
Rcf. 32. 

The remarkable property of the solution is that the 
NT with poor "metallicity" , g < 1 , eventually becomes 
metallic as the length, ft,, of NT is increased. This is not 
the case if the NT is located parallel to the conducting 
gate at distance D <^ h. In the latter case^^^iM one has 
to replace Ch by \n{D/r). 



III. THOMAS-FERMI EQUATION FOR AN 
ARRAY 



Consider an array of parallel NTs located at points, 
Ri, on the substrate, see Fig. 1. To set the Thomas- 
Fermi equation for a given NT, i, in the array, one has to 
take into account that the external potential, leading to 
the charge separation, contains, in addition to eFz, the 
potentials created by charges induced on all other NTs. 
Then the generaUzed Eq. Q reads 

1 

eFz^- p,{z)+ / dz' Vpj(z')5(z,z';R, -Rj), (6) 
9 Jo 



where the kernel, S, is given by 



S{z,z';R) = 



1 



1 

V/(Z + Z')2 + |^' 



(7) 



It is convenient to rewrite Eq. ^ in the "continu- 
ous" form by introducing the position-dependent density, 
p{z, R), and the local concentration of NTs 

AA(R) = ^ 5 (R - RO . (8) 

i 

In the new notations Eq. ((BJ takes the form 

eFz = — p{z, R) 
9 

+ J dR'AA(R')^ dz' S{z, z'; R - R')p{z', R'). (9) 

The Thomas-Fermi equation in the form ^ is convenient 
for further analysis. This is because, as we will see below, 
the distributions, Pi{z), are almost the same for all i even 
for completely random array. 



IV. REGULAR ARRAY 

Consider a regular array in the form of a square lattice 
with a lattice constant, d. Then the coordinates of R, in 
Eq. ijHJ) are {nd, md} with integer m and n. Obviously, 



for the regular array the induced charge density, p{z), is 
the same for all NTs, so that Eq. @ acquires the form 

eF2=(^i±|^^ p{z) + £dz' Se,t{z,z')p{z'), (10) 
where 

Sext{z, z') — 



E 

{m,n}=£{0;0} 



[z - Z'f + [m? + V?) 



\J{z + z' f + {m? + n?) 



(11) 



Here we have isolated the self-action, m = n = 0, of 
a NT. For a single NT this self-action is described by a 
large logarithm Ch- In the array, however, the interaction 
is screened at distances \z — z'\ > d hy the neighboring 
NTs. To account for this screening, the logarithm, Ch, 
in Eq. ifTU)! is replaced by Cd — \n{d/r) . 

It is apparent that both terms in the sum (|ll|l diverge 
due to the contributions from large m, n. However, the 
divergent parts in both terms cancel each other. Physi- 
cally, this reflects the screening by the image charges (see 
Fig. 1). Replacing the sums over m and n in Eq. fTT|l 
by integrals, we obtain the following expression for the 
kernel Sext{z, z') 

Sext{z, z') = 27r7Vo z + z' -\z- z' 

A-kMo\z'Q{z - z') + zQ{z' ~ z) , (12) 

where <d{x) is the step-function. Note that both steps, 
replacing Ch in Eq. (|10(l by Cd and replacing the sums in 
Eq. (|11|1 by integrals are by no means obvious and require 
justification. This justification is provided in the next 
Section. In the present Section we demonstrate that the 
integral equation (|10|l with the kernel (|12|) can be solved 
analytically. Upon taking the derivative from both sides 
of Eq. lfTU|l . we obtain 



eF ■ 



1 + 2g/:„ 
9 



dp 
dz 



+ 47rA/'o / dz'p{z'). (13) 



It is now easy to see that, within a factor, the first term 
in Eq. \i'6l is the second derivative of the second term. 
Thus, Eq. (|13() can be viewed as a second-order differen- 
tial equation with respect to J^dz' p{z'). The solution of 
this equation, satisfying the condition p{0) = 0, has the 
form 



p{z) — pq sinh (z/a) . 
where po is defined as 

eFga 



Po 



(1 + 2gCd)cosh{h/a)' 



(14) 



(15) 
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and the length, a, is given by 



^^1 l + 2gCd 

2 V TTCjAfo 



(16) 



The above expression for the induced charge density con- 
stitutes the main result of the present paper. It is seen 
from Eq. (|14|l that a plays the role of the penetration 
depth of the external field into the array. In the limit of 
very low density, Wo = < 1/ft,^, Eas. ((njl - l(TB|) re- 
produce the result Eq. ^ for a single NT, as it could be 
expected, since the mutual influence of NTs, separated 
by a distance ^ /i is negligible. 

It also follows from Eqs. l(T^ - (fT?))l that in the limit of 
large A/q, such that a <^ h, the induced charge density 
is concentrated near the NT's tips and falls off' towards 
the substrate exponentially, as exp{— (/i — z)/a}. This 
weak penetration of external field into the array is a con- 
sequence of the collective screening. Indeed, in terms of 
screening properties, the array of a high density can be 
viewed as a homogeneous metallic plate. Our crucial ob- 
servation, however, is that the penetration depth exceeds 
parametrically the lattice constant, d, both for large and 
small values of the interaction parameter, g. 

By virtue of the relation a/d > 1, there are many NTs 
within the penetrations depth. This, in turn, suggests 
that Eqs. (|14|I - H16|I apply to the random array with av- 
erage areal concentration of NTs Ao = d~^. 



V. DERIVATION 

In the previous section the derivation was based on two 
intuitive assumptions. Namely, we have assumed the self- 
action of a NT in the array is described by £d instead of 
Cfi for an isolated NT, and that the sum over m and n in 
Eq. 1)11(1 can be replaced by the integral. To justify these 
assumptions, below we calculate the sum Eq. 1)11(1 more 
accurately. In order to do so, we employ the following 
(obvious) identity. Consider a two-dimensional vector, 
b, with projections bx,by. Then, for arbitrary x, we have 



1 



V^2 + |bP 

dqxdqy 
2tt 



exp 



qlx 



exp(i5x6x + iqyby). 



(17) 



by = 



To use this identity in Eq. ((11(1 . we set bx = nd 
md. Then the summation over n and m can be readily 
performed, yielding the sum of (5-functions, i.e. 




(18) 



where p and / assume all integer values. After that the 
rhs in Eq. acquires the form 



dz'p{z')S{z,z'), 



(19) 



where S{z, z') is given by 
27r 



S 



{z,z') = ^(z-i-z'- -I- 



1 



exp 



(t) ^p^V¥^^ 



■ exp 



(20) 



The first term in l(2l)(l describes the continuous limit and 
comes from p = Z = in Eq. 1(181) . It coincides with 
the kernel, Sext{z^z'), defined by Eq. p2(l . The remain- 
ing sum over p and I recovers the kernel S'o in Eq. 
The easiest way to see this is to replace the sums over 
p and I by corresponding integrals, which would imme- 
diately yield Sq{z^z'). However, such a replacement is 
justified only when the large number of terms contribute 
to the sum. This is the case only when the condition 
\z — z'\ < d is met. For \z — z'\ > d the sum over p and I 
in Eq. I(20() is dominated by the terms p = 0, I = ±1 and 
I — 0, p = ±1, and is exponentially decaying function of 
\z —z'\. This suggests that So{z, z') should be substituted 
into Eq. ((19() . in which the integration over z' should be 
restricted to the interval \z — z'\ < d. Within this in- 
terval, p{z') in the integrand of Eq. 1(19(1 can be replaced 
by p{z). The remaining integral yields 2Cd — 21n(d/r), 
similarly to Eq. with h replaced by d. The product 
2Cdp{z) is nothing but the first term in the rhs of Eq. I^. 

The restriction of the integration interval in Eq. ((191) to 
\z ~ z'\ < d for the part of S{z, z'), coming from the sec- 
ond term in Eq. I(2UI) . is, in fact, a delicate step. Although 

this part decays as exp|— 27r|z — z'|/c?| outside this in- 
terval, the behavior of p(z') outside this interval is also 
exponential, namely, it increases as exp(z'/a). There- 
fore, the restriction of the integration interval in ((19(1 is 
allowed only when the exponent in p{z') is slower, i.e. 
a is > d. However, we know from Eq. 1(16(1 that this is 
indeed the case. 



VI. FLUCTUATIONS OF INDUCED CHARGE 
DENSITY IN A RANDOM ARRAY 

The conclusion drawn in Sect. IV that the charge dis- 
tribution Eq. ((14() applies not only to regular but also 
to a random NT array was based on the relation a > d 
between the penetration depth and the lattice constant. 
Thus, this conclusion pertains only to the "body" z a 
of the distribution. Since p(z) falls off exponentially away 
from the NT tip, it might be expected that in the tail 
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(h — z) ^ a the randomness in NT positions would ter- 
minate the applicability of Eq. (|14|l . To verify this fact, 
one can incorporate the positional disorder into Eq. Q 
perturbatively, i.e. to find the correction to the average 
p{z) linear in the fluctuation of the NT density. Then 
the region of applicability of Eq. H14(l to the random ar- 
ray can be established from the condition that the typical 
disorder-induced correction is smaller than the average. 
This program is carried out below. 

In a random array, fluctuations in the areal concen- 
tration of NTs, SU{R) = 7V(R) - (A/'(R)), lead to the 
fluctuations in the distribution of the induced charge den- 
sity Sp{z,'R.). We linearize Eq. @ and in the first order 
over the fluctuations obtain 



H{5p{z,Ii)} = J^(z,R), 
where the integral operator, Ti., is defined as 
7i:{/(z,R)} =4^AAoa2 [/(z,R) 

dR' [ dz'Sext{z,z';Il-R')f{z\Il'] 



(21) 



47ra2 



(22) 



The rhs in Eq. H21|l is the potential created by the fluc- 
tuation, SAfCR), of the density of NTs with unperturbed 
charge distribution, Eq. H14|l . As follows from Eq. 
this potential is given by 



jr(z, R) = - y dR' SU{R') dz'S{z, z'- R - R')p{z') 



(23) 



The fact that the kernel of the integral operator, H, de- 
pends on the difference (R — R') suggests transforma- 
tion from (5A/'(R) and Sp{R, z) to the Fourier harmonics 
(5A/'(q) and (5p(z, q), where q is the in-plane wave vector. 
Upon the Fourier transform, Eq. (|21|l assumes the form 

AttNoo^ 5p{z,q) + / dz'5p{z',c\) 



q Jo 



the following second-order differential equation with z- 
independent coefficients 



where 7^ is defined as 



1 



.F"(z,q)-g2.F(z,q) 



(26) 



(27) 



Note that the rhs of Eq. (|26|1 can be cast in the following 
simple form 

T"{z, q) - q'^T{z, q) 47r^A;'(q)po sinh(z/a). (28) 



It can be now seen from Eq. (|28|l that Ag(5A/'(q) sinh(z/a) 
is a particular solution of the differential equation H26|l . 
However, to flnd the solution of the original integral equa- 
tion (|24|l . one has to complement the particular solution 
with the solution of the homogeneous equation, i.e. to 
write 



Sp(z,q) 



Xq sinh(7qz) -I- A, sinh(z/a) 5Af{q) 
P{z,q)dU{q), (29) 



and find the constants Xq and by substituting Eq. H29|) 
into Eq. This yields 



Xq = 



2po 



cosh{h/a) + qa smh(h/a) 
7q cosh{-fgh) + gsinh(7g/i) 



(30) 



(31) 



x{exp[-|z-z'|g] -exp[-(z + z')g]} =.F(z,q), (24) 
where the rhs is proportional to (5A/'(q) 
Airapo 



Note that Xq diverges at small q. However, the full so- 
lution Eq. H29|) remains finite in the limit q 0. It also 
satisfies the obvious condition Sp{0, q) = 0. 

Equations H29|l - H31|) allow one to quantify the effect of 
disorder in the NT positions on the distribution of in- 
duced charge. The most interesting case is /i 3> a, when 
this distribution is determined by collective screening in- 
volving many NTs. In this limit Eq. (|29|l can be simpli- 
fied by replacing sinh(/i/a) and cosh(/i/a) by exp(/i/a) 
and introducing zi = {h — z) ^ h. Then h drops out 
from zi-dependent part of P(q) in Eq. (|29(l . and we ob- 
tain 



J'{z,q) 



sinh(gz) exp{—qh} (25) 



q{a'^q^ — 1) 

cosh(ft,/a) + qasa\\\{h/a) — agsinh(z/a) \ 5N{q}. 



The structure of the kernel in Eq. H24|) is similar to that 
in the unperturbed equation Eq. (|10|l . It appears that, 
due to this similarity, Eq. H24|l can be solved analytically 
in the same way as Eq. (|10|) . Namely, upon taking the 
second derivative from both sides, Eq. H24|l reduces to 



h/a 



Afoa^ 



(32) 



exp(-7gZi) - exp(-zi/a) 



exp(-7g^i) 
{q + -fq)hq + l/a) 



The form (|32|l is very convenient to study the effect of 
disorder in the "tail", i.e. at large zi. Indeed, assuming 
Gaussian fluctuations in 5Af{R), so that 



(W(qi)5A^(q2)) = 27rAAo5(qi - qa). 



(33) 
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the variance of random fluctuations in the induced charge 
density can be expressed as follows 



{5p[z,f 



i j dR(<5p(zi,R) 
^ |rfqP(zi,q)2. 



(34) 



Here A is the normalization area. It is now seen that the 
g-dependence of P{zi , q) is dominated by the first term in 
Eq. H32|l . The reason for this is the following. As was ex- 
plained in the beginning of this Section, the applicability 
of Eq. l(T^ , obtained for the regular array, is expected to 
be terminated in the random array at "depths" zi that 
are > a. At these depths the average field is strongly 
suppressed. On the other hand, for zi > a one can use 
the expansion 7^ — ^Jq^ + \/a? — \ + ^ • This, in turn, 
suggests that characteristic values of the wave vector, 9, 
are q < 1/(001)"'^/^ < l/o- Then the typical ratio of the 
second and the first terms in (|32|l is q^a? <C 1. It also fol- 
lows from the expansion of 7^ that the main exponents in 
5p{zi) and in the average p{zi) are the same. Upon ne- 
glecting the second term in Eq. (|22J), the q-dependence of 
P{zi,q) acquires the form P(zi, (7) c>c exp(— ag^zi/2)/g^. 
Then the g-integration in Eq. (|34|l can be easily per- 
formed. We will present the final result as the ratio of 
variance, {[6p{h — zi)]'^) and the square of average charge 
density 



[Spjh - z,W 




(35) 



The above result offers the quantitative answer to the 
question about fiuctuations of the induced charge den- 
sity due to the randomness in the NT positions. In 
particular, it can be concluded from Eq. (|35|l that the 
disorder-induced fluctuations in the charge density are 
negligible, if zi < Aoa'^. Since this value is much bigger 
than a, Eq. H35|l conflrms our earlier claim that Eq. H14() 
applies not only for regular, but also for the random ar- 
ray. However, this applicability is limited by the distance 
zi < Aoa'^. For larger zi the variance exceeds the aver- 
age suggesting that the charge density strongly fluctuates 
within the plane zi = const. Note, however, that these 
fluctuations are smooth with characteristic scale (2:10)^/^, 
which is much smaller than zi , but much bigger than the 
penetration depth, a. 

As a final remark of this Section, we point out that 
the lower is the density of the random array the bigger 
is the depth, zi, down to which Eq. H14I) applies, as it 
follows from Eq. (|35|l . However, the magnitude of the 
decay of the charge density, p{h — zi)/ p{h), is governed 
by the ratio zi/a. For zi = Aoa"^, this ratio depends on 
the density of the array only weakly (logarithmically) . 



VII. IMPLICATIONS FOR FIELD EMISSION 

A. Single NT 

It is commonly accepted that the field emission cur- 
rent, J, from the NT tip is described by the Fowler- 
Nordheim law2^ 



HJ{F)/Jo)\ - 



4V2mI^3 



(36) 



where Jq is the prefactor, m is the electron mass, W is the 
work function, which, in principle, is dependent on the tip 
geometrjiSSiSLS. Parameter /3 is the field enhancement 
factor. Various applications of the field emission from 
NTs are based on the fact that /3 is large as a result of 
the NT geometry, more specifically, due to the large ratio 
h/r. The expression for the enhancement factor routinely 
used in the fitting the experimental I-V curves^ is /3 = 
Ch/r, where C ~ 1 depends on specific geometry of the 
tip. Within the Thomas-Fermi description of the induced 
charge distribution, outlined in Sect. II, the expression 
for field at a distance, zi, from the NT tip is given by the 
derivative of the potential, (/'(zi), created by the induced 
charges 

F.ndi^i) = = f dz p{z)S^{z, zi), (37) 

dzi azi Jq 

where p(z) is given by Eq. jSJ. Then the evaluation of 
the integral l|T7l) yields 



Find{zi) 



^Chr 



i{l,r/zi}. 



(38) 



where we had assumed Fimi 3> F. It is seen from Eq. H38|) 
that the enhancement factor indeed has the conventional 
form, /? = Ch/r, with C « {2Ch)^^ for zi < r, but 
it falls off with increasing zi. This suggests that for 
low enough applied flelds, when the electron tunnelling 
length ~ W/Find exceeds r, the I-V characteristics devi- 
ates from the Fowler-Nordheim law. In order to estimate 
this deviation, we substitute 



{zi) = ^ \ -e(r - zi) -f 1 - ln(r/zi) e(zi - r) 
2Lh r 



(39) 



into the tunnelling action 



ln(J(^^)/Jo) - 



dzi ^/W - e(j){zi), 



(40) 



where Zt is the turning point at which the expression 
under the square root is zero. In (|4()|l we had neglected 
the bare potential eFzi . It is now convenient to measure 
the electric fleld in terms of Fq , defined as Fq = W j e(3r = 



7 



2WCh/eh. The integral in Eq. H40() can be reduced to the 
error function, erf {x), after which Eq. (|40|l acquires the 
form 



ln(J(F)/Jo) 



4V2mW^3 



GiFo/F), 



3eh (3Fo 

where the dimensionless function G(r) is defined as 



(41) 



G(t) = t, r < 1; 



1/2 



3 /¥ 



4 V r 



exp(r- 1) erf (Vr- 1), r > 1. (42) 



The plot of the function G(r) is shown in Fig. 2. Strictly 
speaking, the Fowler-Nordhcim region, corresponds to 
T < 1, where the slope of G(t) is identically unity. How- 
ever, G(t) can be linearized around r > 1, where the 
slope is larger. For example, at r = 2 the slope is ~ 2. 
This can be interpreted as a two-times reduction of the 
enhancement factor, /3, in Eq. (|41|l . A significant reduc- 
tion of the enhancement factor (e.g., 30 times) occurs 
around r w 5. It should be noted, that, since Eqs. ||1T)) 
and were derived neglecting the bare potential, their 
applicability is limited by r < Tmax, where r^ax cor- 
responds to the applied field F = Fo/rmax, for which 
the turning point, Zt, in the tunnelling action Eq. H4l)|l 
reaches the value W/F. The latter condition can be 
rewritten in the form {r^nax - 1) = ln(/iTmaa/?'), yielding 
Tmax ~ ln(/i/r) = Ch. For r > T^ax, i-e. for applied 
fields F > Fo/Tj-aax^ the I-V characteristics is given by 
the Fowler-Nordheim law H36() with (3 — 1. 

Overall, Fig. 2 indicates that, for low enough applied 
fields there are significant deviations from the Fowler- 
Nordheim law in the I-V characteristics of an individual 
NT. For such fields the linearity of the Fowler-Nordheim 
plots shows only within very narrow interval of F. On the 
experimental side, there are reports, e.g., Ref. !4d| where 
applicability of the Fowler-Nordheim law was demon- 
strated within a rather wide (exceeding 3 times) interval 
of change of F. In other reported measurements, see, 
e.g., Ref. ^Ij, linearity of In J vs. 1/F holds only within 
a limited (less than 2 times) interval of applied fields. 
Whether or not Eq. H41|) . derived for a single NT, is suit- 
able to fit experimental results depends crucially on the 
geometry of the array, as we discuss below. 



B. Array of NTs 

As it was mentioned in the Introduction, increasing 
the density of NTs in the array leads to dramatic sup- 
pression of the enhancement factor. To illustrate this 
point, consider first the array of low density, when the 
tunnelling length is much smaller than the distance be- 
tween the neighboring NTs. Then the field created by 



induced charges near the tip of a given NT can be calcu- 
lated from Eq. ||2ZI), with p{z) given by Eqs. d)), l(T3)l. 
This calculation yields the generalization of the field en- 
hancement factor Eq. (|38|l to the case of the array of low 
density 



Find{zi) 

F 



tanh(/i/a) 



i{l,r/2i}. (43) 



The above expression recovers the enhancement factor for 
a single NT in the limit a — > oo, or equivalently, A/o — > 0. 
For a < ft., we conclude that the enhancement factor for 
the array, compared to the single NT is suppressed as 



mo) 
m 



hCd 



r2 



1/2 



< 1. 



(44) 



We now turn to the high-density array. In such an array, 
the tunnelling length of an emitted electron can exceed 
the distance, A/q , between the neighboring NTs. Then 
the form of the tunnelling barrier is not anymore given 
by the potential created by a single NT with charge dis- 
tribution modified by neighboring NTs. Instead, one has 
to use a general expression 



(I)a{zi) ^ / dz posmh[{h- z)/a]S{z,zi), (45) 
Jo 

where S is the kernel defined by Eq. H2U|I . The first 
term in the kernel H2U|I corresponds to the "continuous" 
limit, and yields a contribution Fzi to (f>Aizi). The sec- 
ond term in Eq. (|20|l exhibits different behavior for large 

(|z - zil > JVo^^^^ and small (|z - zi| < A^g"^/^) dis- 
tances. The expression for (/)^(zi) that captures both 
long, and short-distance behaviors has the form 



<I>a(zi) 




a 0(zi — d) 



2'Ka -f d 



exp(— 27rzi/(i) — 1 



-e(d- zi)ln ( — ) exp(-zi/a) J.. (46) 



The long-distance behavior is described by the first term 
in the square brackets. It represents the correction to 
the "continuous" first term, Fzi, in (|46|l due to discrete- 
ness of the array. Clearly, the field enhancement due 
to this term is negligible. It is the second term in the 
square brackets that is responsible for the field enhance- 
ment. The physics, captured by this term, is that at 
distance zi < Mq the tunnelling electron "sees" not 
only the NT, from which it was emitted, but also neigh- 
boring NTs. This term, however, does not contain the 
NT radius, r, which was set zero in the derivation of 
Eq. (|46|) . The dependence on r can be reinstated in the 



8 



same way as for a single NT in Eq. H39|l . namely 



e(zi-r)e(rf-zi) 



1 - In I — 

zi 



exp(-[^i-r]/a) , (47) 



where we have retained only the part, responsible for the 
field enhancement. The subsequent calculation of I-V 
characteristics using Eq. H47|l is completely identical to 
the case of an isolated NT. The result can be presented 
in the form similar to Eqs. (|^ . 



ln(J(F)/Jo) 



3eh P{Afo)F, 



■Ga{Fa/F), (48) 



1-1 



1x3/2 



Ga{t) - r 

1 \ ( r } 

—= / du (1 + Inu^r) cxp-^ (u-r — l)}- 

VtJi [ '-a i 

(1 + lnu)exp|--(u- 1)| 



(49) 



where Ut, which is related to the turning point, , as 
Ut = Zr/r, satisfies the following equation 



(1 + Inur) exp| {ur ^ !)}• 



(50) 



The boundary value of the external field, Fa, corresponds 
to the turning point z^- = {Fa/ F)r — tat > r. The field 
Fa is related to the single NT boundary field, Fq, by 



Fa = ^F,. 



(51) 



Therefore, the variable, ta, defined above, is related to 
the corresponding single NT variable, t, as 



TA 



2tt CdAfoh'^ 
1? 



-I 1/2 



(52) 



The function Ga{ta) is plotted in Fig. 2 together with 
the function G{t). We see that, while the parameter 
Moh? changes within a wide interval, the function G{ta) 
remains close to G(t). This means that the form of the 
I-V characteristics for the array is similar to that for a 
single NT. The difference essentially amounts to rescal- 
ing of the characteristic field, Fq, by the factor Cdh/ Cha- 
in other words, increasing the density of the array results 
in suppression of the field emission without the change 
of the shape of the I-V characteristics. It should be 
pointed out that this conclusion pertains to the array 
of randomly positioned, but completely identical NTs. 
As we will see below, the situation changes dramatically 
when the heights of NTs are random. 
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FIG. 2: Thin solid line is the Fowler-Nordheim law 
ln(J/Jo)| oc r, where r = Fq/F. Thick solid line is 
the dimensionless current-voltage characteristics | ln( J/Jo)| vs 
T = Fo/F of an individual NT plotted from Eq. Dashed 
curves (1), (2), and (3) are the current- voltage characteristics, 
In J| vs r = Fa/F plotted from Eq. ^ for (h/r) = and 
dimensionless array densities Moh? = 100, 10, 1, respectively. 



VIII. DISCUSSION AND CONCLUDING 
REMARKS 

The main result of the present paper is Eq. (|14|l that 
describes the crossover of the induced charge density dis- 
tribution from a single NT to the dense regular array 
of NTs. We have also demonstrated that Eq. H14|l ap- 
plies to the random array. Disorder in the NT positions 
terminates the applicability of Eq. H14II only at large dis- 
tances from the tips, where the charge density had al- 
ready dropped significantly. Concerning the field emis- 
sion, our calculations quantify a strong suppression of 
the emission current in a dense array. The field enhance- 
mcnt factor falls off with the NT density, A/q , as Nq 
[sec Eq. (|44|) ]. This conclusion might seem to contradict 
the majority of experiments, where high enhancement 
factors for dense arrays were reported. More precisely, in 
majority of experiments the dependence of emission cur- 
rent on the NT density is simply not addressed, and the 
I-V characteristics are interpreted basing on the proper- 
ties (such as work-function) of individual NTs. In fact, 
in those few papers where this issue is addressed, the 
suppression of field emission with increasing the NT den- 
sity is pointed out on qualitative level. The resolution of 
this contradiction, in our opinion, lies in the fact that in 
realistic situations the heights of NTs in the array are 
widely dispersed. To get an insight how this dispersion in 
heights affects the field emission, consider a regular ar- 
ray in which one NT is higher than others by /ii, which 
is much larger than the average NT separation, Ag , 
but much smaller than h. Within the interval Q < z < h 
the distribution of charge in this "sticking out" NT is 
"enforced" by the neighbors, and is given by Eqs. H14|l - 
(|16|l . However, within the interval h < z < {h -\- Zi) 
this NT "sees" the rest of the array as an equipotential 
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plane. From this observation we immediately conclude 
that, within the interval h < z < [h + zi), the charge 
density in the sticking out NT is given by Eq. © with z 
replaced by (z — h). This, in turn, suggests that the en- 
hancement factor of external field in the sticking out NT 
is high and is equal to hi/2£hir, as follows from Eq. 

The above reasoning suggests that the conjecture that 
the field emission current is dominated by sparse stick- 
ing out NTs allows to account for the high values of 
the enhancement factor observed in experiment. We 
will now demonstrate that this conjecture also allows 
to explain why the dependence ln(J(F)/Jo) follows the 
Fowler-Nordheim law H36|l within a wide interval of F, 
while Eqs. (|41|l and H42(l predict strong deviations from 
Eq. H3t)|) as F is decreased. 

Obviously, the probability, P(hi), to find within the 
array a NT that sticks out by hi decreases with hi. Con- 
tribution of NTs with given hi to the emission current is 
determined by the product 

Jh,(xe^^{[-A{2mWr^f'^/m] G{Fq / F)^P[hi), (53) 

where the first term is the tunneling action, which de- 
pends on hi through Fq = 2W Chi / ehi, the function G is 
defined by Eq. H42I) . Since the tunneling action increases 
rapidly with hi, the product H53|l has a sharp maximum 
at a certain optimal hi. Therefore, ln(J(F)/Jo) is deter- 
mined by the logarithm of the rhs of Eq. (|53l) taken at 
optimal hi. The natural choice for P{hi) is the Poisson 
distribution, exp{—hi/H). We can also use the fact that 
within the interval 2 < r < 10 the function G{t) can be 
approximated with high accuracy by the power law 

G(r) « 0.23r9/2. (54) 

Using this approximation, the optimal hi can be easily 



found analytically. It is convenient to cast the final result 
for \n{J{F)/Jo) in the following form 

where the "effective enhancement factor is defined as 

Firstly, we see from Eq. H55|l that the I-V character- 
istics is very close to the Fowler-Nordhcim law, since 
the exponent 9/11 is close to 1. This should be con- 
trasted to the I-V characteristics of a single NT, for 
which I \n{J{F)/Jo) \ cx F'^/^ as follows from Eq. 
Secondly, the effective enhancement factor (156(1 is large 
and depends rather weakly on the work function W. 
Summarizing, the dense array of NTs can exhibit the 
Fowler-Nordheim field emission provided there is a suffi- 
cient spread in the NT heights. In fact, this conclusion 
is in accord with reported experimental findings. In par- 
ticular, direct imaging of emission intensity by means of 
scanning^i*^ and electron emission^'^-'*'* microscopy re- 
veals that only a tiny portion of NTs in (10^'^ or even 
smaller) contributes to the net current. 
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